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Is gluonic color-spin locked phase stable? 
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We study the gluonic color-spin locked (GCSL) phase in dense two-flavor quark matter. In this 
phase, the color and spatial rotational symmetries are spontaneously broken down to SO(2)di ag 
with the generator being an appropriate linear combination of the color and rotational ones. The 
Meissner masses of gluons and the mass of the radial mode of the diquark field in the GCSL phase are 
calculated and it is shown that this phase is free from the chromomagnetic and Sarma instabilities 
in the whole parameter region where it exists. The GCSL phase describes an anisotropic color and 
electromagnetic superconducting medium. Because most of the initial symmetries in this phase are 
spontaneously broken, its dynamics is very rich. 
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It is plausible that a color superconducting phase is 
realized in the interior of compact stellar objects 
Matter inside compact stars should be in /3-equilibrium 
and be electrically and color neutral. Owing to these 
conditions and the non-negligible strange quark mass, a 
mismatch Sfi between the Fermi momenta of the pair- 
ing quarks is inevitably induced. This is crucial for the 
quark-pairing dynamics Q. 

Moreover, when the mismatch Sfi increases, the Meiss- 
ner masses of gluons become imaginary both in the 
gapped (2SC) and gapless (g2SC) two-flavor color super- 
conducting phases [3(. A similar phenomenon has been 
found also in the three-flavor gapless color-flavor locked 
(gCFL) phase $MM- 

As was shown in Ref. 0, the chromomagnetic insta- 
bility at least in the 2SC/g2SC phase is closely con- 
nected with the appearance of tachyonic plasmons in 
the physical gluonic channels. This supports the sce- 
nario with gluon condensates (gluonic phase) proposed 
in Refs. HQ. 

Another problem in the g'2SC phase is the existence 
of the Sarma instability, which corresponds to the nega- 
tive mass squared of the physical diquark excitation (the 
diquark Higgs mode) at zero momentum. It was also 
found that the diquark Higgs mode has a negative veloc- 
ity squared in the g2SC region [10| . In Refs. [TJ [HJ a 
similar instability is discussed. 

One of the urgent problems in this field is to resolve 
these problems. Besides the gluonic phase @,@|, a num- 
ber of other candidates for the genuine ground state have 
been proposed [H Q EE ES EE M M M M M H3 • 

In Ref. [24) [ . a numerical analysis for two gluonic 
phases, the minimal cylindrical gluonic phase II and the 
gluonic color-spin locked (GCSL) one [9|, l24| . was per- 



formed 1 . It is shown that the gluonic phases are ac- 
tually realized in wide regions of the parameter space 
and they are energetically more favorable than the nor- 
mal, 2SC/g2SC, and the single plane wave Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) [H, [H, HI phases. 

A formula for the Meissner screening mass was devel- 
oped in Ref. [3(| and also the Meissner masses in the 
minimal cylindrical gluonic phase II were examined. It 
was found, however, that the chromomagnetic instabil- 
ity is only partially resolved in the minimal cylindrical 
gluonic phase II. 

In this paper, the stability of the GCSL phase is stud- 
ied. We calculate the Meissner masses for gluons and 
also the mass of the diquark Higgs field at zero momen- 
tum. It is shown that the GCSL phase resolves both of 
the chromomagnetic and Sarma instabilities in the whole 
region where the phase exists. 

As in Refs. [1, HU, in the analysis, the gauged 
Nambu-Jona-Lasinio (NJL) model with two light quarks 
will be used. We neglect the current quark masses and 
the ("i/^) 2 -interaction channel. The Lagrangian density 
is then given by 
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with 



(2) 

where e and e a are the totally antisymmetric tensors in 
the flavor and color spaces, respectively. We also intro- 
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For the earlier works, see Refs. I25l |26H . The extension to the 
model with nonzero temperature was studied in Ref. l27l . 
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FIG. 1: The dynamical solutions of the GCSL phase with 
a s = 1. The values A = 653.3 MeV and fi = 400 MeV were 
used. 



FIG. 2: The mass of the diquark Higgs field in the unit of 
Ml (= 4/i 2 / 7! " 2 )- The values A = 653.3 MeV, = 400 MeV 
and a a — 1 were used. 



duced gluon fields A^, the QCD coupling constant g, the 
generators T a of SU (3) and the structure constants f abc . 
The quark field ip is a flavor doublet and color triplet. 
The charge-conjugate spinor is defined by if) = C^) T 
with C — Z7 2 7°. We do not introduce the photon field. 
On the other hand, the whole theory contains free elec- 
trons, although we do not show them explicitly in Eq. |T]) . 
In /3-equilibrium, the chemical potential matrix /x for up 
and down quarks is fio = /il — /i e Qemj with 1 = l c ® If, 
and Q cm = l c ® diag(2/3, — 1/3)/, where /i and /i e are 
the quark and electron chemical potentials, respectively. 
(The baryon chemical potential \xb is given by /ib = 3/z.) 
The subscripts c and / mean that the corresponding ma- 
trices act on the color and flavor spaces, respectively. 

In the fermion one-loop approximation, the bare effec- 
tive potential including both gluon and diquark conden- 
sates is given by 



with the counter term, 
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where A and S g 1 denote the diquark gap and the fermion 
propagator inverse with gluon condensates in the Nambu- 
Gor'kov space, respectively. We also added the free elec- 
tron contribution. Since the bare potential has a diver- 
gence, a counter term is required. The dimensional regu- 
larization is nice for gauge theories, but not for the NJL- 
type models. We here take into account only differences 
of the free energies with and without the chemical poten- 
tials. Such effects should be physical. This is a similar 
idea for the subtraction scheme considered in Ref. ||. 
Let us define the renormalizcd effective potential by 



V c , = -- 



d 4 P 
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(5) 

In this prescription, even if we use the regularization 
scheme with the sharp three-momentum cutoff A for the 
loop integral, we can remove artificial mass terms of glu- 
ons like A 2 ^* 2 . 

In this paper, we study the GCSL phase with 
an ansatz for the gluon condensates, 

M8 ^ 3 (A«), K = g{Ai)=g{At). (6) 

The dynamics of the GCSL phase was analyzed in 
Ref. [2J|. In order to find the dynamical solutions, we 
searched for minima of the (neutral) effective potential 
imposed the neutrality conditions. We also converted 
the four-diquark coupling constant Ga to the 2SC gap 
parameter A defined at S/j, = and varied the values 
of Ao from the weak coupling regime (Aq ~ 60 MeV) to 
the strong coupling one (Ao ~ f60 MeV). The dynami- 
cal solutions of the GCSL phase are shown in Fig[TJ In 
the analysis, we took realistic values fi = 400MeV and 
A = 653.3MeV. For the gluonic phase, it is required to 
specify the value of a s [= (? 2 /(47r)], because of the exis- 
tence of the tree gluon potential term. In Ref. [3l|, the 
confinement scale is estimated to be not so large. In 
FigfTJ we took a s = 1 as a typical value. We then find 
that the GCSL phase exists in the region 
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76 MeV < A < 161 MeV, 



(7) 
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for a s — 1. While the values of K tend to be small but 
nonzero near the endpoint of the GCSL phase, A ~ 161 
MeV, the gap A and the gluon condensate K are still big 
even in the vicinity of the starting point, Ao ~ 76 MeV. 
(See FigfTJ) These suggest that the first order phase tran- 
sition occurs both at the starting and end points of the 
GCSL phase. 

Note that the values of Ao at the starting point of the 
GCSL phase are rather sensitive to the choice of a s , be- 
cause the tree gluon potential term is non- negligible, i.e., 
as the values of a B are bigger, the windows are wider [24| . 

For reasonable parameter regions, say, 300 MeV < ji < 



500 MeV and 653.3 MeV < A < 1 GeV, the behaviors of 
the solutions are qualitatively unchanged. When we vary 
the values of fj, from 300 MeV to 400 MeV and from 400 
MeV to 500 MeV, the windows where the GCSL phase is 
energetically stabler than the minimal cylindrical gluonic 
phase II are about 15% wider, respectively. On the other 
hand, when the value of A is taken to A = 1 GeV from 
A = 653.3 MeV, the window is about 30% wider. 

We now analyze the Meissner screening masses and 
the mass of the diquark Higgs field in the GCSL phase 
by using the formulae derived in Ref. (30j . 
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where El 2 ... denote the energy eigenvalues of S g 1 . We 
also defined the tree contribution 
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and the transformed vertices 
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and 











(10) 



(11) 



(12) 



in the Nambu-Gor'kov space. For the Dirac's 5-function, 
a lot of expressions are known. In our calculation, we 
use S(x) = lim y/n/irexp(— nx 2 ). We also crosscheck 

n — >oo 

our calculation by using the numerical derivative of the 
effective potential. 



Let us introduce the notations 
d 2 V* 



Mi 



dA 2 
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and 



(M 2 ) ai -b j: n, (M 2 ) ai -b j -c k ,n 



(13) 



(14) 



for the eigenvalues of the Meissner mass (2x2 and 3 x 
3) matrices, where (M 2 ) ai -b^i < (M 2 ) ai -t lj ,2 and so on. 
Notice that there appear mixing terms among the space 
components of gluons owing to the symmetry breaking 
structure of the GCSL phase 0] 



SU(3) C x U{1) 



x SO(3) rot ^ SO(2) diag , 



(15) 



where the unbroken S l O(2) c ii ag symmetry consists of an 
appropriate linear combination of the initial color and 
rotational symmetries. At least the mixing terms which 
exist in the tree contribution should be taken into ac- 
count. To reduce our labor, we might ignore the other 
types of the mixing terms. 
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FIG. 3: The Meissner masses for the GCSL phase without the 
mixing terms in the unit of Mg[= 4a s /i 2 /(37r)]. The values 
A = 653.3 MeV, fi = 400 MeV and a s = 1 were used. 



FIG. 4: The Meissner masses for the GCSL phase with the 
mixing terms in the unit of Mg[= 4a s /i 2 /(37r)]. The values 
A = 653.3 MeV, ju = 400 MeV and a s = 1 were used. 



We depict the numerical results in Figs.|2HI|m the units 
of M| o = 4 ( u 2 /tt 2 and M g = 4a s /i 2 /(3?r). In the whole 
region where the GCSL phase exists, it turns out that 
the diquark Higgs mass and all of the eigenvalues of the 
Meissner masses are positive. We also find that the re- 
lations (M 2 ) lxlx ~ (Af 2 )i^ s (M 2 ) 3x3x s (M*) a . a „ 
(M 2 ) lyly ~ (M 2 ) 3A , (M 2 ) 1-3 . ^ -(M 2 ) lz3 «, 
(M 2 ) 1A ~ (M 2 )^ ~ (M 2 ) 3A ~ -(M 2 ) 3A , 
{M 2 ) 2A ~ (M 2 ) 2A , (M 2 ) 4 ^ ~ (M 2 ) 6 , 6 „ 
(M 2 ) 4 , 4 , ~ (M 2 ) 6x6x ~ (M 2 ) iA , (M 2 ) 5y5y ~ 
(M 2 ) 7b7h , (M 2 ) 5 ,5, =s (M 2 ) 7x7x ~ -(M 2 ) 5s7x and 
(M 2 )g x 8 :c — (-^ 2 )8 i 8 i numerically hold. The behavior 
of (M 2 )8 H 8 B m FigJ3] looks peculiar. Although the nu- 
merical calculation itself is robust, because the behav- 
ior is unchanged by varying the values of q s , /i and 
A, the reason is unclear at present. It is noticeable 
that there are six vanishing masses within the numerical 
precision, {M 2 ) iyAy , (M 2 ) 5x5x , (M 2 ) 6y6y , (Af 2 ) 7 ^ and 
(M 2 ) 4z - 6x! i, (Af 2 ) 5z _ 7x! i. Three among the six should 
correspond to the gauge fixing fields. One might guess 
the other three should be connected with the Nambu- 
Goldstone (NG) bosons owing to the global symmetry 
breaking U(l) x 50(3) ->■ SO(2). (Note that we did not 
introduce the photon field.) This is unclear, however, 
because the kinetic terms are quite important for the ab- 
normal number of the NG bosons 32]. This question will 
be answered elsewhere. 

Since the tree contributions are proportional to K 2 , the 
a s -dependence of the Meissner masses are numerically 
non-negligible. Nevertheless, the results are qualitatively 
unchanged for 0.75 < a s < 1.25. We have also checked 
positivity of the Meissner masses and the diquark Higgs 



mass in the reasonable parameter regions, 300 MeV < 
H < 500 MeV and 653.3 MeV < A < 1 GeV. 

In summary, we analyzed the Meissner screening 
masses and the mass of the diquark Higgs field in the 
GCSL phase. We showed that unlike the minimal cylin- 
drical gluonic phase II the GCSL phase resolves both 
of the chromomagnetic and Sarma instabilities in the 
whole parameter region where the phase exists. In this 
sense, the GCSL phase describes a stable vacuum against 
the fluctuations of the diquark and gluonic channels. 
These results are qualitatively unchanged for realistic 
values of a S) say, 0.75 < a s < 1.25, and for the rea- 
sonable parameter regions, 300 MeV < n < 500 MeV 
and 653.3 MeV < A < 1 GeV. 

Why can the GCSL phase resolve the chromomagnetic 
instability? In the minimal cylindrical gluonic phase II 
only with B = (A 6 Z ) ^ 0, the squared Meissner mass of 
the transverse mode of the 4th gluon becomes negative 
at a certain value of Ao [HJ. Therefore, the GCSL phase 
with (A®) 7^ and (A~.) ^ can be energetically stabler 
than the minimal cylindrical gluonic phase II and resolve 
this instability. Actually, we have found that the GCSL 
phase is energetically more favorable over the minimal 
cylindrical gluonic phase II and the sin gle plane wave 
LOFF phase in a wide parameter region [24|. After the 
rearrangement of the vacuum, the curvatures of the effec- 
tive potential, i.e., the Meissner masses for A\ and A®, 
are likely to be positive. Moreover, unlike in the minimal 
cylindrical gluonic phase II, all of the gluonic fluctuations 
except for A\ and A\ have big tree terms proportional 
to K 2 . Thus the squared Meissner masses tend to be 
positive. 

We also comment on the discrepancy between the 
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starting points of the GCSL phase and the chromomag- 
netic instability of in the minimal cylindrical gluonic 
phase II [13, HO] . The structures of the tree gluon poten- 
tial between the two phases are different. In addition, the 
values of the dynamical solutions B and K are the same 
order, but numerically disagree [24]. It is thus possible 
that the discrepancy numerically occurs. 

The GCSL phase has other noticeable features; (1) 
Both of the chromoelectric and chromomagnetic fields 
are spontaneously generated. (2) The GCSL phase de- 
scribes an anisotropic medium. (3) The medium is a color 
and electromagnetic superconductor. 

In the confinement picture, the symmetry breaking 
structure of the GCSL phase and the cylindrical gluonic 
phase I is the same [9] . How about the relation between 
them? In addition, the question whether or not the ab- 



normal number of the NG bosons occurs is still open. 
Last but not least, is the GCSL phase the global vacuum 
in dense two-flavor quark matter? We hope to return to 
these problems elsewhere. 
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